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Abstract 

There is an ongoing effort to quantify entanglement of quantum pure states 
for systems with more than two subsystems. We consider three approaches 
to this problem for three-qubit states: choosing a basis which puts the state 
into a standard form, enumerating "local invariants," and using operational 
quantities such as the number of maximally entangled states which can be 
distilled. In this paper we evaluate a particular standard form, the Schmidt 
form, which is a generalization of the Schmidt decomposition for bipartite 
pure states. We show how the coefficients in this case can be parametrized 
in terms of five physically meaningful local invariants; we use this form to 
prove the efficacy of a particular distillation technique for GHZ triplets; and 
we relate the yield of GHZs to classes of states with unusual entanglement 
properties, showing that these states represent extremes of distillability as 
functions of two local invariants. 
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I. INTRODUCTION 



The importance of quantum entanglement, both as a resource for quantum information 
processing and as a ubiquitous feature of quantum systems, has become increasingly appar- 
ent over the last few years fiJ-Q. Recent developments in quantum information theory, in 
particular, have stimulated interest in the quantification and manipulation of entanglement. 

For bipartite pure states an essentially complete theory of entanglement now exists , 
though the situation for mixed states is less definite 0. All descriptions of bipartite pure 
state entanglement start with the Schmidt decomposition. It is possible to find orthonormal 
bases {1*)^} and for systems A and B such that we can write the joint state of the 

system in the form 

\*AB)='52Vpi\i) A ®\i) B , Pi>0, £>i = l- C 1 ) 

i i 

These Schmidt coefficients {pi} are uniquely defined by the state \^?ab), and are equal to the 
eigenvalues of the reduced density matrix p A (or equivalently, of p B )] the bases {(i)^} and 
are eigenbases of p A and p B , respectively, so the local density matrices are diagonal 
in this choice of bases. This choice of bases also minimizes the number of terms needed to 
represent \^ AB )- 

For tripartite or multipartite states, there is no equivalent to the Schmidt decomposition 
(0) 0)- Three main approaches to parametrizing tripartite or multipartite entanglement 
have been followed so far. First, one may choose the local bases to put the joint state into 
a standard form. Often these standard forms are intended to generalize some aspect of 
the Schmidt decomposition in the bipartite case [p|-|T5|]. Second, one may try to identify a 



complete set of locally invariant quantities, functions of the state which are invariant under 



local unitary transformations lp| |12l , |i6l -^ , and which uniquely characterize equivalent states. 
The coefficients of a standard form are obviously such quantities, but they may not have 
readily meaningful physical interpretations. Third, one may identify operational quantities, 
such as the number of Greenberger-Horne-Zeilinger (GHZ) triplets or EPR pairs that can 



be distilled from the state by some procedure |TT| , pTipO| . 

In section II we consider a number of proposals for standard forms of three-qubit pure 
states, concentrating especially on those which generalize some aspect of the bipartite 
Schmidt decomposition: the minimal form fl3-15|, the two-term form |T3| , |22|j2T[| , and the 



Schmidt form. This form was given briefly in JET]] and independently in ]T2] . In section III 
we examine it in greater detail. We give an explicit parametrization of the coefficients in 
terms of five locally invariant quantities, and discuss their physical significance. 

We make use of the Schmidt form to prove analytically the reliability of a proposed 
distillation technique for GHZ triplets from general three-qubit pure states |TT[ ; we present 
this proof in section IV. 

In H Linden and Popescu proposed characterizing the entanglement properties of three- 
qubit states by examining the "orbits" of the states under general local unitary transforma- 
tion. This was carried a step further by Carteret and Sudbery ||19|| , who proved that most 



states have a certain generic behavior under such transformations, but identified classes of 
'special' states which they speculated to have unusual entanglement properties. In section 
V we briefly review these 'special' states, then analytically evaluate the yield of GHZs under 
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the distillation protocol of JTT| and section IV. By examining the yield of these states as 
a function of the invariant parameters from section III, and also of the locally invariant 
"residual tangle" tabc of Coffman, Kundu and Wootters |H| , we verify that these classes of 
states are indeed exceptional by this operational measure, representing extremes of distilla- 
bility or undistillability. We briefly compare the results using this protocol to the recently 



discovered optimal distillation method of [plf| . 
Our conclusions are summarized in section VI. 



and find that they are entirely consistent. 



II. REVIEW OF STANDARD FORMS FOR THREE-QUBIT STATES 

Two qubits can always be represented in their Schmidt decomposition (ffl) 



\^) = y /p\O0) + y Jl-p\ll), (2) 

characterized by a single Schmidt coefficient p (or equivalently 1 — p). Without loss of 
generality, we adopt the convention that p > 1/2 and the corresponding eigenvector is |0). 

Most attempts to define a standard form for a three-qubit state attempt either to gen- 
eralize some property of the bipartite Schmidt decomposition, or make use of the Schmidt 
decomposition between one of the bits and the other two, or both. For instance, we can 
make a Schmidt decomposition between qubit A and qubits B and C, writing the three-qubit 
state in the form 

= VpMaMbc + ^p\l) A \ipi) BC . (3) 

Choosing the Schmidt basis for qubit A guarantees that the correlated states of qubits B 
and C must be orthogonal: (V'olV'i) = 0. 

The sixteen real parameters to describe a generic pure state of three qubits can be reduced 
to fifteen by normalization, and to five which are invariant under the ten-dimensional group 
of local unitary transformations fl8[,|19| ; unfortunately, no single choice of five quantities has 
proven completely satisfactory. One simple parametrization that has been proposed [j§HT0ll 
is the Linden- Popescu-Schlienz (LPS) standard form. One begins with a state in form @. 
One can then choose one of the two correlated states, say \i^o) BC , and find its corresponding 
Schmidt bases. The resulting state for the three qubits has the form 



|V>> = VpWa «|oo) bc + Vi^lii) 



BC 



+ V / T^|1)^ ( 7 (v / T^|00) i3C - a|ll) BC ) + /|01) BC + <7|10) BC ) , (4) 

where p, a and / are real positive numbers, g is complex, and 7 = (1 — f 2 — \g\ 2 Y^ 2 - Together 
these give five independent real parameters. 

The vectors ^0)23 an d I ^1)23 span a two-dimensional subspace of the Hilbert space for 
qubits 2 and 3. It's possible to make an interesting variation on the LPS idea using a result 
of Niu and Griffiths, who showed [2~I] that any such two-dimensional subspace can be given 
basis vectors of the form 



i Xo ) = ^ioo> 2 3+ v /r^|ii> 23 , 

|Xi) = v / ^|01) 23 + v / r^ : |10) 2 3, (5) 
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for some choice of a product basis for the 4-D Hilbert space of the two bits, where q and r 
are real numbers between and 1. Using this basis leads to a unique standard form 



y/p\0) (ay/q\00) + aJl-q\ll) + by/f\01) + bVT=r\10)) (6) 



+y/l-p\l) (^-6*v^|00) + av/f|01) + a v / r^7|10)J , 

where a is real and a 2 + \b\ 2 = 1. This then gives five independent real parameters: p, q, r, 
a, and the phase of b. This form treats the |0) and |1) terms more symmetrically than LPS; 
however, there is still a lack of symmetry under interchange of the bits. 

More interesting from a fundamental point of view are attempts to generalize some 
aspect of the Schmidt decomposition. Three such properties suggest themselves. First, the 
Schmidt decomposition is the choice of orthonormal bases for the local Hilbert spaces which 
minimizes the number of terms needed to represent the state. Second, any two qubit state 
can be written as the sum of only two product vectors. (For N- dimensional systems, N 
product vectors are needed.) Third, the Schmidt decomposition diagonalizes the reduced 
density matrices of the local subsystems. No single representation for tripartite systems has 
all three properties, but they can be generalized individually. 

Acfn et al. |13j have shown that all three-qubit states can be written in the form 

= A |000) + Aie^|100) + A 2 |101) + A 3 |110) + A 4 |lll) (7) 

by a suitable choice of basis, where the Aj are all real and positive and is a phase between 
and 7T. With only five terms, this is a minimal description, and in that sense a generalization 
of the bipartite Schmidt decomposition. A similar form has been described by Higuchi and 
Sudbery |T4T |, 



= A e^|000) + Ai|100) + A 2 |010) + A 3 |001) + A 4 |lll) (8) 
which has the added benefit of being symmetric under interchange of the qubits. Carteret, 



Higuchi and Sudbery |15j have shown how to generalize this construction to give a unique 
minimal representation for systems of any dimension. These minimal forms have practical 
benefits: with a small number of terms, they can simplify the calculation of locally invariant 
quantities f23| . However, the Aj and <j) themselves have no obvious physical interpretation. 

This minimal property can be generalized in another way, by relaxing the requirement 
that the product vectors be orthogonal. Acin et al. and Diir, Vidal and Cirac have also 
shown [ ]T3"| , p2| that almost all three qubit states can be written in the form 

\^abc) = Hilaibicx) + yu 2 e l</i |a 2 & 2 c 2 ), (9) 

where the vectors are normalized but not orthogonal. There are six real parameters, fi 2 , 
(ai\a 2 ), (bi\b 2 ), (ci\c 2 ) and </>; imposing normalization reduces this to five. 

Interestingly, not all three qubit states can be written in the form (PI); a small subclass 



of states require a minimum of three product terms [13]. Diir, Vidal and Cirac made use 
of this result to prove that there are two classes of three-qubit pure states which cannot 
be interconverted with nonzero probability p2[. The class that requires three terms is a 



three-parameter family, and is characterized by vanishing residual tangle tabc — PI (see 
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section V). Acm, Diir and Vidal also used this form to demonstrate a method of converting 
a single copy of a three qubit state into a GHZ triplet with maximum probability . 



The third generalization is to find bases for all three qubits which diagonalize their 
reduced density matrices. That is, one can simultaneously put each bit in its Schmidt de- 



composition with respect to the other two. This form was proposed in [JTTJ] and independently 
The state has the form 



m 



12 



a|000) + 6|001) + c|010) + d|011) 
+e|100) + /|101) + <?|110) + fe|lll), 



(10) 



which looks just like a generic three-qubit state with 16 parameters. However, using each 
of the three qubits in turn we can write in a form similar to (^), with orthogonality 
conditions which impose restrictions on the possible values of the coefficients in (flOf) . We 
can use these relationships to reduce these coefficients to five independent parameters, as 
we show in the next section. 



III. PARAMETRIZING THE SCHMIDT FORM 

By redefining the relative phases of the basis vectors 

lO^liy,- exp(%)|0} ) exp(^)|l), (11) 

we can choose to make four of the coefficients real. A convenient choice is to make a, d, f, g 
real, while b, c, e, h remain complex. The state must also be normalized, which imposes the 
condition 

a 2 + \b\ 2 + \c\ 2 + d 2 + \e 2 \ + f + g 2 + \h\ 2 = 1. (12) 

This leaves 11 undetermined parameters. 

We can now express the larger eigenvalues Pa,b,c of the reduced density matrices Pa,b,c 
in terms of the coefficients: 

PA = a 2 + \b\ 2 + \c\ 2 + d 2 , etc., (13) 

(the smaller eigenvalues obviously being 1 — Pa,b,c)- Finally, the states \if>o,i) kl correlated 
with basis vectors |0) ■ and |1) . must be orthogonal to each other. This gives three more 
equations: 

ae* + bf + cg + dh* = 0, etc. (14) 

Because these equations are complex, they are equivalent to six real equations. 

Combining these restrictions, we now have fourteen equations in sixteen unknowns. Thus, 
in addition to the eigenvalues Pa,b,c we would expect there to be two more free parameters. 
Can we identify reasonable candidates for these parameters? It turns out that natural choices 
are the two probabilities a 2 and \h\ 2 . These parameters are symmetric under interchanges of 
the three qubits, and have a fairly simple physical interpretation: they are the probabilities 
of all three qubits giving the same result (0 or 1, respectively) when measured in their 
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Schmidt bases. Moreover, the coefficients of the other state vectors can all be calculated in 
terms of the five probabilities a 2 , |/i| 2 , and Pa,b,c, u P to a sign. 

Define p sum = Pa + Pb + Pc- The expressions for the norms of the coefficients are then 
simple: 

2 2 {2p C ,B,A ~ ^)\h\ 2 ~ (Psurn ~ PC,B,A ~ l)(2a 2 - p sum + 1) 



2psum 3 



2 2 (%PA,B,C ~ l)^ 2 - feum ~ Pa,b,c ~ 1 ) ( 2 1 /z | 2 +p s 



d,f,9= ^ ^ • (15) 

The phases of b, c, e are more complicated. If we define the variables <pb,c,e by b = 
\b\ exp(z</>f,), c = |c|exp(i0 c ), and e = |e|exp(i0 e ), the constraint equations ( JUty imply after 
a bit of algebra that 

cos(0 6 ), sin(0 6 ) = {Q 1;2 /\b\){T2adf + g{a 2 + d 2 + f 2 - g 2 )), 
cos(0 c ), sin(0 c ) = (Q lt2 /\c\)(T2adg + f(a 2 + d 2 - f 2 + g 2 )), 
cos(0 e ), sin(0 e ) = (Q 1>2 /\e\)(T2afg + d(a 2 - d 2 + f + g 2 )), 
cos(0 ft ), sin(0 h ) = (Q lj2 /\h\)(T2dfg + a(-a 2 + d 2 + f 2 + g 2 )), (16) 

where Qi and Q 2 are two constants. We can solve for the values of Qi and Q 2 by using the 
identity sin 2 (0) + cos 2 (0) = 1 and substituting flT5| ) for |6|, . . . ,g. 

In the Schmidt form for three-qubit pure states, each of the five parameters has a reason- 
ably straightforward physical interpretation. The three parameters Pa,Pb,Pc are the larger 
(i.e., p > 1/2) eigenvalues of the reduced density operators for each of the three qubits, and 
correspond to the probabilities of obtaining the more likely of the two possible outcomes 
(which by convention we label |0)) when we measure each of the qubits in its Schmidt basis. 



These parameters are closely related to the minimum absolutely selective information (IT 
for each qubit, which is given by the entropy function 

min Si = -(pi\og 2 pi + (1 - p 4 ) log 2 (l - p^). (17) 

This quantity is the minimum amount of fundamentally unpredictable classical information 
generated by carrying out a measurement on qubit i, given a free choice of measurement 
basis fllTf . By using the Schmidt form to choose measurement bases we can simultaneously 
minimize the absolutely selective information for all three qubits. 

The parameters pa,Pb,Pc range from 1/2 to 1 (since they are defined to be the larger 
eigenvalues of their corresponding local density matrices). Similarly, a 2 ranges from to 
1, and \h\ 2 from to 1/2. However, this does not mean that these parameters can take 
arbitrary values within these ranges. Some choices of parameter values correspond to no 



physical state, and give nonsensical values for (]T5|) and (jig). 

In particular, the local probabilities must obey the triangle inequalities 

p A (l - p A ) +p B (l- Pb) > Pc(l - Pc), 
PbO- - Pb) + Pc(l - Pc) > - Pa), 

p c {l - Pc) + Pa(1 - Pa) > Pb(1 - Pb); (18) 

these imply, for instance, that if Pa — 1 then ps = Pc- The restrictions on a 2 and \h\ 2 
are more complicated, but they too display an interdependency in their range of values. In 
particular, as p S um ^ 3 we must have a 2 — > 1 and \h\ 2 — > 0. 
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IV. PROOF OF DISTILL ABILITY 



The Schmidt form can provide analytical insight when addressing specific problems. 
For example, the efficacy of a recently proposed tripartite distillation protocol [JTTJ can be 
demonstrated with its help. 

Consider a state of three qubits in an arbitrary product basis, which can be written in 
the form fllQl) . We can straightforwardly calculate the quantity p A (l ~ Pa) 

p A (l — Pa) = \ a f — be\ 2 + \ag — ce| 2 + \ah — de\ 2 

+ \bg-cf\ 2 + \bh-df\ 2 + \ch-dg\ 2 } (19) 

This expression is a polynomial in the coefficients and their complex conjugates, and is 
correct in any basis. If the state is in the Schmidt form, this simplifies to 

p A (l - PA ) = (a 2 + \b\ 2 + \c\ 2 + d 2 )(\e\ 2 + f + g 2 + \h\ 2 ). (20) 

Let us assume that we have written the state in Schmidt form, such that the states 
{|0), |1)} for each qubit j are eigenstates of the local density matrix with eigenvalues pj and 
1 — Pj, respectively. Suppose we now perform a weak measurement on each of the three 
qubits. First, allow each qubit to interact with a separate ancilla bit initially in state |0), 
such that 

|0> ® |0) anc - y/T=l\0) ® |0) anc + v^|0) ® |l) anc , 
|l)®|0) anc -|l)®|0) anc , (21) 

where e<l. Then measure the three ancilla bits. With a probability of e(p sum ) one will 
find one or more of the ancilla bits in state |l) anc , in which case the procedure has failed. 
Otherwise, this step has succeeded and the three qubits are now in a new state with slightly 
different coefficients a', b', . . . , h'. The changes in the coefficients are 

Aa = -(e/2)(3 - p Bnm )a, 
A{b, c, e) = -(e/2)(2 - p sum )(b, c, e), 
/, g) = — (e/2)(l - p mm )(d, f, g), 
Ah = (e/2)p snm h. (22) 

This very simple form results because the state is in Schmidt form. After this procedure the 
bases for the three bits will generally no longer be the correct Schmidt basis (though it will 
be close to it), so the expression (p0|) cannot be used; but (|19"D is always correct. Thus we 
get a change m pa(1 — Pa) 

A[p A (l-p A )} = -e(A-2{p A + p B +p c )){\af -be\ 2 + \ag-ce\ 2 ) (23) 
-e(3 - 2(p A + p B + Pc ))(\ah - de\ 2 + \bg - cf\ 2 ) 
-e(2 - 2(p A + p B +Pc))(\bh - df\ 2 + \ch - dg\ 2 ) 
= -e(3 - 2(p A +p B + Pc))pa(1 - Pa) 

-(e/2)(\af - be\ 2 + \ag - ce\ 2 - \bh - df\ 2 - \ch - dg\ 2 ) 
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By making use of equations (p~3|) and (0), this expression simplifies to 



A\p A (l- PA )} 



which using ([ 

A\p A (l- PA )] 



(2(p A +p B + Pc) - 3)pa(1 - Pa) 
+p A (a 2 - \e\ 2 + \h\ 2 - d 2 ) + d 2 - a 2 

further simplifies to 

e(2p A - 1) 



(24) 



t 2(a 2 + \h 2 \)(p B +p c -V 
2p A + 2p B + 2p c - 3 L 

-{2pa - l)(p A +p B + Pc - 1)(pa + Pb + Pc - 2) 



(25) 



The prefactor to (^) is strictly positive, as is the first term inside the brackets. The second 
term is positive if p A + p B +Pc < 2; any state that satisfies this criterion will evolve towards 
the GHZ state and have a nonzero yield. 

For p A + pb +Pc > 2, the sign of ( Effi) depends on the relative sizes of the first and second 
terms inside the brackets. The last two equations of (p~5| ) show that for p A + p B + pc > 2, 
the fact that f 2 + g 2 > implies 

2a 2 {p B +p c -l) > (2p A +p B +pc-2){p A +p B +p c -2), (26) 

which yields the inequalities 

2a 2 (p B +pc-l)- (2p A - l)(p A + p B + pc - l)(p A + PB +pc ~ 2) 
> {2p A + p B +p c - 2)(p A + p B + p c - 2) 

-(2p A - l)(p A + p B + pc - 1)(pa + Pb + Pc ~ 2) 
= (1 - p A )(p A +p B +p c - 2)(2p A + 2p B + 2p c - 3) > 0. 

This straightforwardly implies 

A[p A (l-p A )} >e(2p A -l)(l-p A )(p A +p B +p c -2) >0. 

Because of the symmetry of the protocol, p B (l—p B ) andpc(l — Pc) must also increase. So one 
step of this protocol must move the state towards the GHZ with nonvanishing probability, 
and will (in general) produce a nonzero yield of GHZ triplets. 

There are three circumstances in which this result can fail. First, no product state can 
ever be distilled to a GHZ by this method. At least one of P A ,Pb,Pc must equal 1 in this 
case, which causes the rate fl25D corresponding to it to vanish. This is not immediately 
obvious from the form of (p5|) , but it is easily checked using (|15D and (|18D — if p A = 1, then 
Pb = Pc — a? \ an d (p5|) is equal to zero. 

Second, there are states with p A + p B + Pc = 2 for which a 2 = \h\ 2 = 0, again making 
(f2"5p vanish. These are a subset of the triple states discussed in section V below, which are 
equivalent to states of the form (|31~D; these states have vanishing residual tangle. Finally, 
it is possible for a state with p A + p B + pc > 2 to evolve to one of these triple states. 
All such states will also have vanishing residual tangle [ 2"2"] , and conversely all states with 
vanishing residual tangle will evolve under this distillation protocol to a triple state with 
Pa + Pb + Pc — 2, and hence have zero yield of GHZs. This can be clearly seen in Fig. 1. 



(27) 



(28) 
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V. ENTANGLEMENT AND DISTILL ABILITY 



Linden and Popescu M proposed characterizing three-qubit states by the dimensions of 
their orbits under the action of the local unitary group. Generically, tripartite pure states 
of qubits have ten-dimensional orbits, equal to the dimension of the local unitary group. 
The very interesting results of Carteret and Sudbery give a complete classification of 



all states for three qubits which behave nongenerically under local unitary transformations; 
these 'special' states have stabilizers of nonzero dimension, and hence orbits of dimension 
< 10 (see [[TjJ). This behavior suggests that these 'special' classes have unusual entanglement 
properties, which might be evident in other measures of entanglement. 

We have numerically simulated the distillation of generic states by the algorithm de- 
scribed above in section IV, in order to determine the yield of GHZ triplets as a function 
of various parameters, especially the parameters used to describe the Schmidt form. We 
have also calculated analytical expressions for the yield of states in the exceptional classes 
enumerated by Carteret and Sudbery. We find that these states are indeed exceptional by 
this operational criterion, as we describe below. 

Most important in calculating the yield of GHZs is the sum of the local eigenvalues 
p snm = Pa + Pb + Pc- This quantity determines the probability of failure in one step of the 
infinitesimal distillation procedure of section IV, with the probability of failure being ep sum . 
If it takes N steps to become sufficiently close to a GHZ triplet, the expected yield is 

Y = ]]_(!- ep sum (n)) w exp |-^ep sum (n)| , (29) 

n=l { n J 

where p sum (n) is the value of p sum at the nth step. In the limit of infinitesimal steps the 
sum inside the exponent becomes an integral. Calculating p sum (n) analytically is no simple 
matter for a general state; the equations (|2~2"D for the change in the coefficients become 
differential equations in the limit, but must be supplemented by an additional change of 
basis between steps, since in general the bases will no longer be the Schmidt basis for the 
new state. While this is simple to do numerically, analytically it is challenging. 

Fortunately, the classes of exceptional states are generally expressible in simple forms 
which make it possible to integrate the equations ( p2"D in closed form, and derive simple 
expressions for the yield of GHZs. Interestingly, the steps of the GHZ distillation technique 
commute with local unitary transformations. Because of this, the distillation procedure 
preserves the stabilizer of the initial state, and hence must take 'special' states to other 
'special' states of the same type. This gives another way of understanding why certain 
special states are not distillable to GHZs. 

In addition to a 2 , \h\ 2 , and p snm , we looked at the dependence of the GHZ yield on one 



other locally invariant quantity. This is the residual tangle of Coffman et al. which can 
be written 



r, BC = 2(AfAf + AfAf), (30) 



where A*/ and \% are the (positive) eigenvalues of the matrix \jpijpij- Here p^ is the 
density operator for the two-party ij system, and is the "spin-flipped" density operator: 
Pij = {p y ® a y)P*iji. a y ® a y)- ^ has been suggested [JTJ] that the residual tangle is a measure 
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of the irreducible three-way ("GHZ-type") entanglement of a tripartite state, beyond any 
two-party ("EPR-type") entanglement that may be contained in such a state. As such, it is 
of particular interest in discussing distillability below. Also, its square tabc 2 is a polynomial 
quantity, which makes it analytically tractable. 

Triple States. For this set of states the residual tangle vanishes flSfl. We previously 



described states in this set as "triple" states [|IT|], because they are equivalent under local 



unitary transformations to states with just three components: 

|Vtr) = 6|001} + c|010) + e|100). (31) 

Carteret and Sudbery |1| refer to these as "beechnut" states; they all have p sum > 2. For 
triple states with p sum > 2 each step of the infinitesimal distillation protocol reduces p sum , 
but leaves the state a triple state. If p sum = 2, the actions on the three qubits cancel 
out, leaving the state unchanged. States of this type have vanishing primary yield for the 



tripartite distillation protocols described in section IV and in indeed, Diir, Vidal and 
Cirac have shown that no procedure can transform one copy of a state with zero residual 
tangle into a GHZ with nonzero probability p2 . 



Because the distillation procedures of section IV and |TT[ preserve the classes of 'special' 
states, it is easy to see why they cannot produce GHZs from triple states; because all triple 
states have p sum > 2, they cannot include the GHZ state (p sum = 3/2) as a limit. The 
set of product states (or "bystander states" in the terminology of Carteret and Sudbery) is 
similarly undistillable. The result of Diir, Vidal and Cirac, however, goes beyond this, since 
it assumes nothing about the symmetry of the procedure. 

The symmetric version of state ( |3~1~D (with b = c = e = 1/V3) is termed by Diir, Vidal 
and Cirac the "W" state, and seems to fill a role for the zero residual tangle states similar 
to the role filled by the GHZ for all other states: it is, in some sense, maximally entangled. 
We will say a bit more about this below. All other 'special' classes include the GHZ as a 
limit, and therefore are distillable. 

Generalized GHZ states. These states can be written in Schmidt form 

\ip) = a|000) (32) 

They have pa = Pb = Pc — fl2 > residual tangle tabc — 4a 2 /i 2 . A single step of the 
infinitesimal distillation procedure gives a new generalized GHZ with coefficients a' = a+Aa, 
ti = h + Ah: 

Aa = -(e/2)(3-jw)a, Ah = +(e/2)p sum h, (33) 
so (|29|) can readily be evaluated to give the yield of GHZs 



Y = 1 - VI - tabc = (2/3) (3 - Psum ) . (34) 

These states are the most distillable three qubit states as a function of both tabc and p S um; 
we can see this in Figures 1 and 2 below. 
Slice states. In Schmidt form these are 

|^) = a |000) + d|011) -e|100) ae = alh, (35) 
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plus similar states derived by permuting the order of the bits. These states have Pb = Pc = 
a 2 + e 2 , pa = a 2 + d 2 , tabc — 4(a/i + de) 2 = 4a 2 (h + e 2 /h) 2 . Imposing normalization and 
the orthogonality condition on fl35|) we see that this is a two-parameter family of states. 
For these two parameters we may choose a 2 and h 2 , or equivalently pa and ps- One step 
of the infinitesimal distillation protocol applied to state (|35| ) leaves qubits B and C in their 
Schmidt bases, but not qubit A; a change of basis must be applied to A to put the new state 
in Schmidt form. This new state is still a slice state, and has new parameters p' A = Pa + ^-Pa, 
Pb=Pb + Ap B , 

Ap A = 2e(p A + 2p B - l)p A , 

Ap B = 2e(p A + 2p B - 2)p B - ep A (p A + Pb - l)/(2p A - 1). (36) 

The yield is difficult to evaluate analytically, but numerical evidence shows that generic 
slice states are not extremes of distillability. With each step of the distillation protocol, the 
parameters Pb = Pc approach 1/2, but pa actually moves away. However, when p B = Pc = 
1/2, this subclass of slice states does have extremal behavior. Carteret and Sudbery term 
this subclass the maximal slice states. 

Maximal slice or Slice-ridge states are of form (E5|) with a 2 + e 2 = 1/2 = ps = pc- This 
subclass is parametrized by a single number, which can be taken to be p A - Because only pa 
is larger than 1/2, there is no need to perform the GHZ distillation procedure on qubits B 
and C; performing it on A alone preserves the form of the state, with 

Ap A = -(e/2)p A (l-p A ), (37) 

giving a yield of GHZs 



Y = 1- VI ~ tabc = 2(2 -p snm ). (38) 

The expression for the primary yield in terms of the residual tangle is identical to that for 
the GHZ-type states, while in terms of p sum it is not. In terms of tabc it is one of the most 
distillable types of state (see Fig. 1). In terms of p sum (Fig. 2) it appears to be one of the 
least distillable types of states; this is because maximal slice states have the minimum tabc 
of all states with a given p sum . 

In addition to these 'special' states, there are two classes of states that deserve additional 



attention. While these states have stabilizers of zero dimension p5J like generic states, these 
classes are also preserved by the above distillation protocols. Like the 'special' states, they 
extremize distillability SIS db function of tabc and p sum . 

Generalized triple or Tetrahedral states. These states can be written 

\if>) = 6|001) + c|010) + e|100) + h\lll). (39) 

We are mainly interested here in the symmetric state b = c = e; for this case Pa = Pb = 
Pc = 2b 2 , Psum — 2. This form is preserved by the steps of the infinitesimal distillation 
protocol, which make the coefficients evolve according to (22); the yield is easily integrated 
according to (p9|) to give Y = 2(2 — p &ma ) = 4(1 — 3b 2 ); the residual tangle is tabc = 
166 3 \/1 — 3b 2 = \J(A — Y) 3 Y/27. This yield is identical to that of the maximal slice states 
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as a function of p sum , but not as a function of tabc] from Figs. 1 and 2, we see that they 
are states of minimal distillability in terms of both p snm and tabc- 

Zero residual tangle (ZRT) states. Diir, Vidal and Cirac have shown that no states with 
tabc — can be converted to GHZ triplets with nonzero probability, so Y = 0. They also 
showed that all such states can be written in the form 

|^) = a |000) + 6|001) + c|010) + e|100). (40) 

This is in general not in the Schmidt form of section III. These states include the triple 
states a = as a subclass (for which ([40]) is in Schmidt form). The triple states form a 
boundary of this set, and any ZRT state will evolve under the distillation protocol to a triple 
state. These states have p sum > 2. 

All these 'special' states have symmetries which account for both their enlarged stabilizers 
and their extremal distillability. One way of seeing this is to note that the various standard 
forms given in section II, which for generic states all require distinct bases, often coincide 
for these special states. For instance, the generalized GHZ states (^) are simultaneously 
in Schmidt, two-term, minimal, LPS and Griffiths-Niu form. ZRT and triple states cannot 
be written in two-term form, but can be written with three terms; the triple states ([H|) are 
simultaneously in Schmidt, minimal, three-term and Griffiths-Niu form. Slice states written 
in the form fl35f) are simultaneously in both Schmidt and LPS standard forms. 

It is easiest to see how the the distillability of these states compares to that of generic 
states by plotting their yields Y as a function of p sum and tabc along with the numerical 
results for a large sample of randomly generated states. We have plotted these quantities in 
Figs. 1 and 2, with the families of 'special' states indicated. We see that most of these states 
are indeed special as far as distillation is concerned: they form the boundaries of the plotted 
regions. The quantity tabc does seem to be closely related to distillability, as conjectured, 
though this relationship is not exact; for a given value of tabc states with a range of Y values 
exist, but the range is not very wide. This range is bounded at the top by the generalized 
GHZ and maximal slice states, and at the bottom by the symmetric generalized triple state. 
All ZRT states have Y = tabc — 0. There is also a relationship between p sum and Y, 
though again for a given p sum there is a range of Y values. This range too is bounded above 
by the generalized GHZs, and below by the ZRT states, generalized triples and maximal 
slice states. These upper and lower bounds are both linear; the upper bound is exactly the 
same as that for Bernstein and Bennett's Procrustean technique ||, reflecting the fact that 
generalized GHZ states can be distilled by exactly the same techniques which work in the 
bipartite case. 

A reasonable question is to what extent these yields are artifacts of the particular distil- 
lation protocol we use. After all, this technique is only one possible way of producing GHZs, 
in general not the optimal method even for a single copy of a three-qubit state. 

Fortunately, we can actually answer this question. Recently, Acm, Jane, Diir and Vidal 
PH have discovered the optimal algorithm for transforming a single copy of a three-qubit 
state into a GHZ. This involves performing a POVM on each of the three bits, designed to 
project the states in the two-term representation onto tri-orthogonal vectors. Finding the 
correct POVM for an arbitrary state involves maximizing a somewhat involved function, 
but is easily done numerically. We have done so for a large sample of random states, as well 
as for the members of the 'special' classes enumerated by Carteret and Sudbery. 
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The optimal yield is higher, in general, than that of the infinitesimal algorithm of section 
IV, though they are surprisingly close for most states. However, for the 'special' states, the 
yields are identical. In other words, the infinitesimal distillation technique gives the optimal 
yield for these classes of states. Quite remarkably, if we plot Figures 1 and 2 for the optimal 
GHZ distillation protocol, the figures look completely unchanged. 

Thus we can see that by both the optimal and the infinitesimal techniques, these classes 
of special states extremize the yield of GHZs as a function of both tabc and p S um- This 
strongly supports the conclusion that these states do indeed have unusual entanglement 
properties, and are worthy of further study. 



VI. CONCLUSIONS 

We have examined tripartite entanglement from both an analytical and an operational 
point of view. In the bipartite case, which is well understood and to which we have turned 
for clues, the analytical and operational aspects of entanglement are closely related: the 
entanglement properties of a single copy are given by the locally invariant parameters, the 
Schmidt coefficients, which also determine their operational characteristics. We have looked 
for similar connections in the three-qubit case. Here at least five locally invariant parameters 
are required, as opposed to just one in the two-qubit case. We have examined several ways 
of choosing these five parameters, looking in particular at generalizations of the bipartite 
Schmidt decomposition. One representation in particular, the "Schmidt form," has useful 
properties which made it simple to prove the efficacy of the infinitesimal GHZ distillation 



protocol of [IT |; it can also be parametrized in terms of five physically meaningful quantities. 

We have looked for connections between these parameters and yields in distilling GHZ 
triplets, as well as connections with the residual tangle of Coffman et al. We have shown that 
the 'special' classes of states enumerated by the theorem of Carteret and Sudbery extremize 
the distillation yield as functions of the residual tangle tabc and p sum = Pa + Pb + Pc- 

Although a certain amount amount of progress towards understanding tripartite entan- 
glement has been made, at least for qubits, many important questions remain unanswered. 
For example, the number of states in the asymptotic minimum reversible entanglement gen- 
erating set (MREGS) p0| , ^6| , p7| for three-qubit states, and for tripartite states in general, 



is still unknown. No asymptotically reversible (or optimal but irreversible) distillation tech- 
nique for GHZ states is known. The search for solutions to these and related problems is 
ongoing. 



ACKNOWLEDGMENTS 

We would like to thank H.A. Carteret, W. Diir, R.B. Griffiths, A. Sudbery and G. Vidal 
for many useful conversations. This work was supported by NSF Grant No. PHY-9900755. 



13 



REFERENCES 



C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres and W.K. Wootters, Phys. 
Rev. Lett. 70, 1895 (1993). 

D. Deutsch, A. Ekert, R. Jozsa, C. Macchiavello, S. Popescu and A. Sanpera, Phys. 
Rev. Lett. 77, 2818 (1996). 

C.H. Bennett, H.J. Bernstein, S. Popescu and B. Schumacher, Phys. Rev. A 53, 2046 
(1996). 

M.A. Nielsen, Phys. Rev. Lett. 83, 436 (1999). 

G. Vidal, Phys. Rev. Lett. 83, 1046 (1999); D. Jonathan and M.B. Plenio, Phys. Rev. 
Lett. 83, 1455 (1999). 

M. Horodecki, P. Horodecki and R. Horodecki, |quant-ph / 960700~9 . 

A.V. Thapliyal, Phys. Rev. A 59, 3336 (1999). 

N. Linden and S. Popescu, Fortsch. Phys. 46, 567 (1998). 

J. Schlienz, Ph.D. thesis. 

O. Cohen, unpublished. 

O. Cohen and T.A. Brun, Phys. Rev. Lett. 84, 5908 (2000). 
A. Sudbery, J. Phys. A 34, 643, (2001). 

A. Acin, A. Andrianov, L. Costa, E. Jane, J.I. Latorre and R. Tarrach, Phys. Rev. Lett. 
85, 1560 (2000). 

A. Higuchi and A. Sudbery, Phys. Lett. A 273, 213 (2000). 

H. A. Carteret, A. Higuchi and A. Sudbery, J. Math. Phys. 41, 7932 (2000). 
N. Linden, S. Popescu, A. Sudbery, Phys. Rev. Lett. 83, 243 (1999). 

J. Kempe, Phys. Rev. A 60, 910 (1999). 

V. Coffman, J. Kundu and W.K. Wootters, Phys. Rev. A 61, 2306 (2000). 
H.A. Carteret and A. Sudbery, J. Phys. A 33, 4981 (2000). 

C.H. Bennett, S. Popescu, D. Rohrlich, J. A. Smolin and A.V. Thapliyal, to appear in 
Phys. Rev. A, |quant-ph/990W3;. 



A. Acin, E. Jane, W. Diir and G. Vidal, Phys. Rev. Lett. 85, 4811 (2000). 
W. Diir, G. Vidal and J.I. Cirac, Phys. Rev. A 62, 062314 (2000) 
A. Acin, A. Andrianov, E. Jane, R. Tarrach, |quant-ph/0009107j . 



Chi-Sheng Niu and Robert B. Griffiths, Phys. Rev. A 60, 2764 (1999). 
H.A. Carteret, private communication. 

N. Linden, S. Popescu, B. Schumacher and M. Westmoreland, |quant-ph/99 12039 . 
S. Wu and Y. Zhang, |quant-ph/ 0003020 . 



14 



Figure 1. Here we plot the primary yield of GHZ triplets from the infinitesimal distillation 
algorithm of section III vs. the square of the residual tangle tabc 2 f° r various 'special states' 
as well as a random sample of generic states. We see that all states lie between two curved 
boundaries; the generalized GHZ and maximal slice states lie on the upper boundary while 
the generalized triple states lie on the lower boundary. The triple states all have both tabc 
and the yield equal to zero. Interestingly the maximal slice states appear to be high-yield 
states when plotted against tabc, but low-yield when plotted against p sum = Pa + Pb + Pc', 
for a given value of tabc these states minimize p sum . 



Figure 2. Here we plot the primary yield of GHZ triplets from the infinitesimal distillation 
algorithm of section III vs. p sum = Pa+Pb+Pc for various 'special states' as well as a random 
sample of generic states. We see that all states lie between two linear boundaries; the 
generalized GHZ states lie on the upper boundary while the maximal slice and generalized 
Triple states lie on the lower boundary and the triple states are the zero-yield states between 
p snm = 2 and p sum = 3. The upper linear boundary corresponds to the yield of Bernstein 
and Bennett's Procrustean method of EPR distillation in the bipartite case. 
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Pa+Pb+Pc 
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